Maximal Invariants Over Symmetric Cones 
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Abstract 

In this paper we consider some hypothesis tests within a family of Wishart distributions, 
where both the sample space and the parameter space are symmetric cones. For such testing 
problems, we first derive the joint density of the ordered eigenvalues of the generalized Wishart 
distribution and propose a test statistic analog to that of classical multivariate statistics for test- 
ing homoscedasticity of covariance matrix. In this generalization of Bartlett's test for equality 
of variances to hypotheses of real, complex, quaternion, Lorentz and octonion types of covari- 
ance structures. 



1 Introduction 

Consider a statistical model consisting of a sample space X and unknown probability measures Pg, 
with 9 in the parameter space 0. In this setting, the statistical inference about the model is often 
concentrated on parameter estimation and hypothesis testing. In the latter, we typically test the 
hypothesis Hq : 6 € 0q c vs. the hypothesis H : 9 € 0\0o. Systematically, this means to find a 
suitable test statistic t from X to a measurable space Y, and the distribution of tPg, the transformed 
measure with respect to Py under t. 

For a Gaussian model, where the sample space is R", the probability measures are multivariate 
normal distribution N„(0, X) and the parameter space is PD„(R), the cone of n x n positive definite 
matrices over R, some classical examples of such hypotheses are the sphericity hypothesis: 
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for some real k x k matrices A, B, with 2k = n; and the quaternion structure hypothesis 
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for some real mxm matrices A, B, C, D, with Am - n. Incidentally, in each of these hypotheses the 
parameter space in the full model is PD n (R), the parameter space in the submodel is a subcone of 
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PD„(R), and the hypothesis is invariant under an action of a subgroup of the orthogonal linear group 
of R'\ 0(?i, R). 

More generally, suppose O is a symmetric cone and 

M = {Pa- 6 P(X) : cr 6 Q.} (1.4) 

is a statistical model with parameter space O, and 

Mo = {/V 6 P(*) : ctgQo} (1-5) 

is a submodel of At, where the parameter set Oo is a (symmetric) subcone of Q. In this set up, we 
wish to test the hypothesis 

H :o-eQ vs. H : cr e Q\Q . (1.6) 

In the presence of a group action, testing the hypothesis (11.61) leads one to the study of a maximal 
invariant with respect to a group lfT3l[T9l [3ll23l. 

The maximal invariants that arise in testing problems like those stated in (11.11 ). (11.21 ) and (11.31) 
are indeed functions of the eigenvalues of the sample covariance matrix (3l[T][T3l[T9l. Therefore, 
in most situations, the problem is reduced to finding the distribution of the ordered eigenvalues of 
a distribution (often Wishart) over the cone of positive definite matrices over the real, complex, or 
quaternion fields. 

In this respect, studies on the distribution of the ordered eigenvalues of a random matrix with 
values in a classical symmetric cone overlap with similar studies in Random Matrix Theory (RMT) 
EE] [TBI- Therefore, applications of such studies are not merely limited to the testing problems we 
mentioned. Other applications of the distribution of the ordered eigenvalues of a random matrix, 
specially a Wishart matrix, can be found in physics, e.g. lfl8l . principle component analysis, e.g. 
OSnH, and signal processing, e.g. ED USES. 

In this paper we discuss the statistical analysis of models of the form (11.41) . There are many mod- 
els that fall in this category, such as the Wishart, inverted Wishart, beta and hyperbolic distributions. 
In our analysis, we extensively exploit the algebraic and geometric structures of symmetric cones, 
described by Jordan algebras and Riemannian symmetric spaces, respectively (see lfl2l l9l [TTll ). In 
the statistical inference of these models two problems of interest are the estimation of the parameter 
cr, and testing a hypothesis of the form (11.61) . Here we focus mainly on the latter. 

In approaching these problems, we begin in Section |2] with giving a short review of symmetric 
cones and their analysis based on ||9). Essential parts of this section are the classification of irre- 
ducible symmetric cones, or equivalently, simple Jordan algebras, and Pierce decomposition. In 
Section [3] we study special functions over symmetric cones. Our main goal in this section is to 
define zonal polynomials over a symmetric cone. Sections HI [5j and [6] present our contribution to 
the theory of maximal invariants over symmetric cones, introducing a new family of non-central 
Wishart distributions over symmetric cones, deriving the joint distribution of the ordered eigenval- 
ues of such distributions, and proposing an analog of Bartlett's test for testing homoscedasticity of, 
roughly speaking, the covariance matrix. 
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2 The structure of symmetric cones 



In this section we provide the reader with a condense review of the structural analysis of symmetric 
cones which is essential in understanding the generalized notions, such as hypergeometric functions, 
zonal polynomials and the Wishart distribution, we discuss later. More thorough presentation of 
these topics can be found in [9] and elsewhere in this volume. 

2.1 Symmetric cones 

Let (J, (; •)) be a Euclidean vector space and let GL(J) be the set of regular linear transformations 
on J. A convex cone O c / is a set satisfying ax + By e Q. for any a,B > and x,y e O. Let O 
denote the closure of Q in J. The convex cone Q. is proper if O n -O = {0}, full if O - Q. = J. A 
homogeneous cone Q. is a full proper open convex cone with this property that the automorphism 
group of Q, defined by G(Q.) = {g € GL(7) : gQ. = Q}, and consequently the connected component 
of the identity in G(Q), denoted by G, act transitively on O. The homogeneous cone O is called 
symmetric if it is homogeneous and self-dual, i.e., Q.* = {y e V : (x\y) > 0,Wx € Cl\{0}} = D.. 
Every symmetric cone is the finite product of irreducible symmetric cones, where a symmetric cone 
is irreducible if it is not the direct product of two or more symmetric cones. A familiar example of 
an irreducible symmetric cone is PD r (R). 

A well-known property of a homogeneous cone Q is that the stabilizer of each x e Q., G x , is a 
maximal compact subgroup of G. For a symmetric cone Q, if we set K = G n 0(7), where 0(7) 
is the group of orthogonal transformation of /, then K is the stabilizer of an element of Q. and, 
consequently, a maximal subgroup of G. 

Another interesting property of a homogeneous cone Q, is that if we define <p(x) - J a , exp(- (x\y))dy 
for each x e Q., then <p{x)dx is a G-invariant measure on O. This follows from the fact that for any 
geGwe have tp(gx) = |det g\~ <p{x). 

2.2 Euclidean Jordan algebras 

A Jordan algebra J over R. is a real vector space equipped with a multiplication satisfying xy = yx 
and x(x 2 y) = x 2 (xy) for any x, y in J. An associative bilinear form B on 7 is a bilinear form 
satisfying 

B(xy, z) - B(x, yz) for any x, y, z e 

A Jordan algebra with no non-trivial ideal is called simple. Typical examples of Jordan algebras are 
so-called special Jordan algebras. If A is an associative algebra, then the special Jordan algebra A + 
is the vector space A equipped with the Jordan product xoy = (xy + yx)/2. In particular M r (R), the 
set of r x r matrices, and S r (R), the set of r x r symmetric matrices, are special Jordan algebras with 
the Jordan product. 

A finite dimensional Jordan algebra J over R is called Euclidean if there is an associative inner 
product on J. Every Euclidean Jordan algebra has a multiplicative identity element e. The quadratic 
representation of J is the map P : J — > End(7), where P : x i-> 2L(x) 2 - L(x 2 ) and L(x) is the left 
multiplication map by x, i.e., L(x){y) - xy. 
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2.3 Spectral decomposition 

A non-zero element c in J is idempotent if c 2 = c. An idempotent is primitive if it cannot be written 
as the sum of two idempotents. Two idempotents c\ and ci are orthogonal if c\C2 = 0. A Jordan 
frame is a maximal set of orthogonal primitive idempotents c\,...,c r . If J is simple, then for any 

Jordan frame c\, . . . , c r we have ci H h c r = e and the number of elements in this Jordan frame, r, 

is invariant for every Jordan frame. This common number is called the rank of J. Moreover, for each 
x in J there is a Jordan frame c\,..,,c r such that i;\C\ + • • • + % r c r = x, where £1, . . . , % r are called 
the eigenvalues of x. The trace and determinant of x, denoted by tr(x) and det(x), respectively, are 
the sum and product of all eigenvalues of x. 

2.4 Euclidian Jordan algebras and symmetric cones 

For a Euclidean Jordan algebra J the set of squared elements x 2 where x e J is a closed cone. The 
interior of this cone, denoted by / + , is a symmetric cone and is called the cone of positive elements 
of /. Equivalently, / + is the set of all x 2 such that x is regular, i.e., det(x) + 0. Conversely, if Q is a 
symmetric cone and e e Q., then one can construct a Euclidean Jordan algebra J such that J + = Q. 
and e is the identity element of /. 

2.5 Classification of irreducible symmetric cones 

The irreducible symmetric cones are in one-to-one correspondence with simple Euclidean algebras, 
which are classified into four families of classical Jordan algebras together with a single exceptional 
Jordan algebra. 

The first three families of classical Jordan algebras are matrix spaces. More specifically, let 
D = R, C or the quaternion H. Denote by x the conjugate of x in D, %x the real part of x, and H r (D) 
the set of r x r Hermitian matrices over D. Recall that X* , the adjoint of a matrix X, is obtained by 
taking the conjugate of each entries and then transposing the matrix. A matrix X is Hermitian if it 
is equal to X*. This space equipped with the Jordan product X o Y = (XY + YX)/2 and the scaler 
product (X\Y) = %tr(XY) is a simple Euclidean Jordan algebra of rank r and corresponds to the ir- 
reducible symmetric cone of positive definite matrices over R, C or H, respectively. The dimension 
of H r (D) over R is r + r(r - l)d/2, where d, the Peirce constant, is equal to the dimension of the 
space D over R. When D is the set of real numbers H r (R) is indeed S r (R) and S r (R)+ is PD r (R). 

The fourth class of Jordan algebras is a Minkowski space. This is the vector space R x R' I_1 , 
n > 2, equipped with the product (£, x)(i;, y) = + x ■ y,£y + %x). It corresponds to the Lorentz 
cone L„ = {(£,x) e R x : £ > ||x||(. 

The exceptional Jordan algebra can be described as follows. First notice that HxH with the 
product (xi,yi)(x2,y2) - (xi*2 ~ yiyi^yv^i + y2*i) is a non-associative algebra, called octonion 
O. Any element in O can be written asi + jy, where j = (0, 1) and x,j e H. The conjugation in 
O is defined by x + jy = x - jy. The exceptional Jordan algebra is H-3(Q) and its cone of positive 
elements is denoted by PD3(0). 

The following table summarizes the information about simple Euclidean Jordan algebras (9j|6]|. 
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Table 1 : Classification of simple Jordan algebras 
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2.6 Additional properties 

In the remainder of this section, and throughout the rest of the paper, we will assume that Q. is an 
irreducible symmetric cone of positive elements of a simple Jordan algebra / of dimension n, rank 
r and the Peirce constant d defined by n - r + r(r - \)d/2. 
a) For each xeQ. and geGwe have 

1. detL(x) -det(x) n/r , 

2. detP(x) - (detx) 2n/r , 

3. detgx = det(g) r/ " det(x) for each g e G, 

4. detPOXx) = (det;y) 2 detx, 

5. {gx)~ l = g*~ l x _1 , where g* is the adjoint of g, 

6. /Xx)- 1 - PCx" 1 ), 

7. P(x)* = P(x), i.e., P(x) is Hermitian. 

^) ^ acts transitively on the set of primitive idempotents and the set of all Jordan frames. 

3 Special functions on symmetric cones 

The rich geometric structure of symmetric cones has naturally motivated many research studies in 
the field of harmonic analysis. In this section we give a short description of special functions on 
symmetric cones. 

3.1 Gamma function 

Let fix a Jordan frame c\,...,c r in J. For each 1 < j < r we define the idempotent ej - c\ + 
■ ■ ■ + Cj. In general, it can be shown that if c is an idempotent element of J, then the only possible 
eigenvalues of the linear transformation L(c) are 0, 1/2, 1. One can see that the eigenspace, J(ej, 1), 
corresponding to the eigenvalue 1 of L{ej) is a Jordan algebra with the multiplication inherited from 
J. Let Q.j be the cone of positive elements and det (;) the determinant with respect to this Jordan 
algebra. Let Pj : J — > J(ej, 1) be the orthogonal projection on J(ej, 1). The principal minor, Ay(x), 
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is a homogeneous polynomial of degree j on / defined by Aj(x) = der-"(Py(jc)). We extend this 
definition as follows. For each s = (ji,...,j r )eC r we set 

A s (x) = A l (xy i - S2 A 2 (xY 2 ~ S3 ■ ■ ■ A r (x) Sr . 

For each s e C the gamma function is Tq(s) = exp{-tr(x)}A v (x) det(x)~^dx. This integral is 
absolutely convergent if % sj > (j - l)d/2, for j - 1 r. Moreover, 

r j 

r n (5) = (27r) 2 ?nr(^-(j-i)-). 

In particular, identifying z e C with (z, . . . , z) e C, we have 

Tn(z) - I exp{-tr(x)}det(x) z ~'<ix 
is absolutely convergent if %z > n/r - 1. If this is the case, then 

r 

Tn(z) = (27T) 2 ? ]~ [ T(z - 0' - l)d/2). 

7=1 



3.2 Beta functions 

The beta function on a symmetric cone O is defined by the integral 

Bo{a,b)- J A a _ a (x) A fo _ s (e - x)<ix, 
nn(e-n) 

where a, ft e C, and <?-Q = {<?-x:xe Q}. This integral converges absolutely if %aj > (j - l)d/2 

and %bj > (j - l)d/2. In this case B^{a, b) - _£L^_£L. _ ; an( j 

T n (a + ft) 

A a _ 5 (»A fo _«(x - y)dy = Ba{a,b)A a+h _!i{x). 

Cln(x-Cl) 



3.3 The space of polynomials 

A partition is a finite sequence of non-negative integers A - (A\ , • • • , A r ) in decreasing order 
A\ > • • • > A r . The weight of A is \A\ - A\ + • • • + A r . If |/l| = k, then we say /I is a partition 
ofk. 



A function / : J — > R is a polynomial on / if there is a basis {vi, . . . , v n } of J and a polynomial 
p e R[*i, ...,t„] such that for any linear combination x - £" =1 &Vj, 

/(*)=/*£,..., 

One can check that this definition is independent of the choice of basis. The set of all polynomials 
over / is denoted by P(J). The action of G on / can be naturally extended to an action on P(J) by 
defining gp(x) - p(g~ l x). Let ¥^(J) be the subspace of P(7) generated by polynomials gA^, g eG. 
Then every p in P^(J) is a homogeneous polynomial of degree 
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3.4 Spherical polynomials 

Recall that K, the stabilizer of the identity e, is a compact Lie subgroup of G, thus there exists a 
Haar measure on K. For each partition A the spherical polynomial CD,; is 



Jk 



A A (kx)dp K (k), 



where fi^ is the normalized Haar measure on K. The function <S> A is indeed a homogeneous poly- 
nomial of degree \A\ and is invariant under the action of K, i.e., <& A {kx) - Q A (x) for any k e K 
and x e J. Moreover, the spherical polynomial (fr A is, up to a constant factor, the only ^-invariant 
polynomial in f A {J). More precisely, if p is a A'-invariant homogeneous polynomial in P A (J), then 

f p{kx)dfi K (k) = p(e)® A (x). (3.1) 
Jk 

Consequently, for any g £ G 

® A (gkx)dfi K (k) = ® A (ge)® A (x). (3.2) 



X 



IK 

If x e Q and ^Ry > (r - l)d/2, then for any y in O and g e G we have 

f e- tr(x),) ® A (gx) detixY'-rdx = T n (^ + y)A^(y)0 A (gy- i ). (3.3) 
Jn 

3.5 Zonal polynomials 

The Pochhammer symbol for Q. is defined by 

r n (s + A) 

(s) A = ———, jeC. 
Tq(5) 

This definition generalizes the classical Pochhammer symbol. The zonal polynomial Z A is a homo- 
geneous ^-invariant polynomial on J defined by 

Z A (x) = d A ^-® A {x), (3.4) 
{j)a 

where d A is the dimension of the vector space P A (J), and <!>,! is the spherical polynomial. Zonal 
polynomials are A'-invariant polynomials normalized by the property 

tr{xf = V Z A (x) xeQ. (3.5) 

\A\=k 

Notice that the function p(x) = tr(x) k on Q, is a ^-invariant homogeneous polynomial in P A (J). It 
follows from Equation (13.21) that for each x e J and y e O 

f Z A (P(yhkx)dfi K (k) = Z ^ Z ^ X \ (3.6) 
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3.6 Hypergeometric functions 

Let a\,...,a p and b\ b q be real numbers with a, - (i - l)d/2 > 0, bj - (j - l)d/2 > and 

x,y e J. The hypergeometric function p F q is defined by 

, , , v v (^^••■Wj^W^W 

p F q (a u ■ ■ ■ ,a p ,b u . . . ,b q ,x,y) = > > — — — ——. 

For y = e we set p F q {ci\, . . . ,a p ,b\, . , . ,b q ,x) - p F q {a\, , . . ,a p ,b\, . . . ,b q , x, e). This means that 

„, , , , v 1 V ( ai ^ " -( a p)AZ A (x) 

p F q (au . ... a,, fc, ....*,,*) = ^ ^ —. 

This series converges absolutely if p < q and diverges if p > q. Furthermore, o^oW = and 
qF^x) = det(e-x)~ b . 

Suppose x e Q and g e G. By using the integral in Equation (13.61 ) for p < q we get 

I p F q (a\, . . . ,a p ,b\, . . . ,b q ,g(kx))dp K (k) - p F q (ai, . . . ,a p ,b\, . . . ,b q ,x, ge). (3.7) 

JK 

Similarly, if %y > (r - l)d/2, y e Q. and p < q, or y e e - Q. and p = q, then by applying Equation 
(1331) we get 



I e tr(xy) p F q (ai, . . . ,a p ,b\, . . . ,b q ,gx,z)d&t{x) y > dx = 

r n (r) det(j)" r / , + i J F (? (ai , . . . , a p , y, b x , . . . , b q , gy~ l , z). 

We record the following Proposition, later used in deriving the beta type distribution, from 191 . 

Proposition 3.1. Suppose that p < q + 1, > (r - l)c?/2 and ^772 > (r - q)d/2. If g e G such 
that geeOn(e-fi), then 

I p F q (ai, . . . ,a p ,b { , . . . ,b q ,gx)det(x)' 1[ ~~ det(e - x) m ~~dx (3.8) 

Jcin(e-n) 

= B n (qi,q 2 ) P +iF q+ i(a l ,.. . ,a p ,r]ub\, ■■■ ,b q ,m + q 2 ,ge). 

Suppose y e J and x e Q. We define x -ky = P(x?)y and Eig(y|x) equal to the r-tuple of the 
ordered eigenvalues of y with respect to x, in decreasing order. Notice that these eigenvalues are the 
roots of the polynomial p{€) = det(£x - y). In particular, for x = e, we have Eig(y) = Eig(y|e) is 
the ordered eigenvalues of y. Note that for a special Jordan algebra, such as the space of symmetric 
matrices or Hermitian ones, the operation x ★ y is x^yxi. Here we list some interesting properties 
of ★ operation. 

Lemma 3.1. Let x and y be in Q. Then 

a) x-ke = e-kx = x, x-k x~ x = x~ l ★ x = e and (x ★ y)~ l = x~ l ★ 

b) x ★ y and y ★ x have the same eigenvalues, therefore, any ^-invariant function on O returns the 
same value on both x ★ y and y ★ x. In particular, det(;t ★ y) = det(y ★ x) = det(x) det(y), tr(x ★ y) = 
tr(y ★ x) = tr(xy) and Z^x ★ y) = Z^(y ★ x). 
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Proof, a) We only prove the last equality. The other equalities are immediate from the definition 
of x ★ y. From property (5) in Subsection 12. 61 we get 

(x-kyT 1 = {P{x^yT l = P{x^)*~ y~ l = P(x~hy~ l = x~ l * y' 1 . 
b) For an r-tuple (£i , . . . , £ r ) with %\ > ■ ■ • £ r > 0, we define 

U : : &) ' -<iv' 

By this notation, one can check that for any xefiwe have Eig(^^ 1 ) = Eig(x)" 1 . Clearly, the ordered 
eigenvalues of x ★ y are equal to Eig(y|x _1 ). On the other hand, Eig(y|x) is equal to Eig(x|y) _1 (if ( 
is an eigenvalue of x with respect to y, then l/€ is an eigenvalue of y with respect to x). Therefore, 
Eig(x +y)= EigCylx" 1 ) = EgOT^) -1 = Eig^" 1 ★ x" 1 )" 1 = Eig((y- 1 * x" 1 )" 1 ) = Eig(v * x). □ 

Remark 3.1. There is worth mentioning an interesting fact about the x-ky operation on Q. Suppose 
jU and v are two K invariant measures on Q. These can be lift up, respectively, to K bi-invariant 
measures fj* and v* on G. By projecting the convolution measure /r * on fl we obtain a measure 
on Q , denoted by p. ★ v and called the convolution of p, and v. Now if x and y are two random 
vectors in Q. with K invariant distributions p and v, then the distribution of the random vector x * y 
is the convolution p * v ifTOll . 

4 The non-central Wishart distribution 

In multivariate statistics, the so-called non-central Wishart distribution is the multivariate analog 
of the non-cenual chi-square distribution and, similarly, it arises out of the sampling distribution 
of a sample statistic, namely the sample covariance matrix of a multivariate normally distributed 
population. More generally, if X is an N x r random matrix such that X ~ Nnxi-(M, In ® 2) with S 6 
PD r (R), then the density of 5 = X T X with respect to the standard Lebesgue measure is proportional 
to 

det(X)"T exp{-itr(S- 1 5)}exp{-itr(A)} F 1 (i^V, ^£L- l S)l POr(m (S). 

Thus the non-central Wishart distribution is the transformed distribution with respect to a multivari- 
ate normal distribution under the quadratic mapping X i-» X T X. We can generalize this situation as 
follows. 

Let J be a simple Euclidean Jordan algebra, and (E, (-De) a Euclidean space. A (Jordan) repre- 
sentation ip of J over E is a linear mapping ip : J — > End(£) such that tA(x 2 ) = <AM 2 f° r eacn x e ^- 
The representation ifj is called self-adjoint if the linear transformation tfr(x) is self-adjoint for every 
x e /. 

Assume ifj is a self-adjoint representation of 7 on the Euclidean vector space E. For each veE, the 
function x i-» (t/Kx)v|v)£ is a linear form on £\ Therefore, there exists a quadratic map Q : £ — > J 
such that <i/<x)v|v>£ = (x\Q(v))j. The fact that <2(v)|x 2 ) 7 = ||t/Kx)v||| > implies that Q(y) e D, 
for each v e E. 

Proposition 4.1. Let v ~ Nip., <A(cr)) be a normally distributed random vector in E with mean 
vector p € E and covariance ifr(cr) e PD(E), where ifr is a self-adjoint representation of the simple 
Jordan algebra / over E, cr e O, and > 2(n - r). Then the density of x = Q(\) with respect to the 
canonical Lebesgue measure on Q is 

p x (x) = (2h n (^-)det(o-)Tr) exp {-l t r(e)} ■ det(x)^-? Fi(^, \{o-- 1 ★ e) (4.1) 
\ 2r / 2 2r 4 
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where e = cr 1 * Q(p). 

Proof. First we compute g(z) = E[exp{-(z, Q(v)) j}], the Laplace transform of Q(\). For each z e / 
we have g(z) is equal to 

((2tt)t det(^(o-))^ 1 j^xp{-(if f (z)v,v) E }exp{-^(o-r i (v -fi),(y - n)) E }dv 

= ((2tt)t det(cr)^)" 1 • J^exp{-i (2<^(z)v,v> £ + (^(o-y l (v - //), (v - n)) E )\dv. 

We rewrite the expression 2{\]/(z)v, v) E + (i]/(cry l (y - //), (v - as 

<iA(2z + o-' l ){v - iff(2z + (r- 1 )- 1 ^- 1 )^,^ - i/f(2z + c^rVCo- -1 )^ 
+ (-A(2z + o-- l )-^(o-- l )ii, Ho-- l )ii) E + (^(o-- l )fi,iu) E , 

and substitute it in the original equations, obtaining 

g(z) = 2"* det(cr)-£ det(z + jo- -1 )^ exp{-itr(e)} F (^((r- 1 ★ e) ★ (z+ j "" 1 )" 1 ), 

where %{z + ^cr~ l ) > 0. The density of x = Q(\) in Equation (14.11) then follows from the inverse 
Laplace transform of g(z) (see H). □ 

In light of Proposition 14.11 Equation (14.11 ) can be taken as the density of the non-central Wishart 
distribution over the irreducible symmetric cone Q. 

Definition 4.1. A random vector x with values in an irreducible symmetric cone Q. is said to have 
the Wishart distribution with parameters rj > n - r, cr e CI and e e Q., if its probability density with 
respect to the standard canonical Lebesgue measure on / is given by 

exp{--tr(e)} F 1 (-77, ^(cr -1 ★ e) * x)wq(x\t], o-)Iq(x), 

where 

wa(x\r], cr) = r exp{- - trier' 1 x)} tet(x)*^ (4.2) 

This is denoted by x ~ , Wq(j}, cr, e), where rj, cr, and e are, respectively, the shape, (multivariate) 
scale, and non-centrality parameters. In particular, if e = 0, the distribution is called the (central) 
Wishart distribution over a symmetric cone 0, denoted by 

x~W n (ri,o-), 

otherwise, it is called the non-central Wishart distribution over a symmetric cone. 
Remark 4.1. If x ~ "Wsiirj, cr), then its probability density is wq,(x\t], cr). 

Remark 4.2. In accordance with classification of irreducible symmetric cones, the non-central 
Wishart distribution over an irreducible symmetric cone can be categorized as : i) the real type 
"W,-(r], cr, e) over PD, (R), ii) the complex type CW r (rj, cr, e) over PD r (C), Hi) the quaternion type 
WW r {r], cr, e) over PD r (H), iv) the Lorentz type UW n (j], cr, e) over the Lorentz cone L„, and v) the 
exceptional type O'WQ], cr, e) over the octonion cone PD3(0). 
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To explore some of the properties of the non-central Wishart distribution, resembling those of 
the classical non-central Wishart, first we compute the characteristic function of the non-central 
Wishart distribution. By definition this is equal to (p x (z) = E[exp{/tr(xz). Thus we have 

<Mz) = (23"T n (|)det(0-)2") 1 exp{-itr(£)} J^expl-itrCCo-- 1 - 2iz)x)} 

■ oFi(-ri, -(a _1 ★ e) ★ x)d&t{x)^ r '~~'dx 



= (25^det(o-)^) _I exp{-itKe)}det( - - 2lZ )~^ 

1 1 1 cr" 1 - 2iz , 1 
• 2 71, 4 ( 2 } * ( * 

Simplifying this we get the characteristic function of x is given by 

i n 1 1 -1 

<f> x (z) = det(e - 2iazY n expj--tr(e)} e\p{-tr(e - 2iaz) e)}. 

Using this we can show that the convolution of two independently distributed Wishart is a Wishart 
distribution again. More precisely we have: 

Corollary 4.1. Suppose x ~ "Wsi(r]i, cr,e\) and y ~ 'Wofrqi, cr, ei) are independently distributed. 
Then x + y ~ "Waim + cr, e\ + e 2 ). 

Proof. The characteristic function of x+y is the product of the characteristic functions of x and y. It 
remains to check that this product is equal to the characteristic function of "W^qx +772, cr, ei+e£). □ 

Other properties of the non-central Wishart distribution on an arbitrary irreducible symmetric 
cone are similar to those of the classical non-central Wishart distribution. Here we list some of these 
properties that can be shown to hold as in Proposition 3.1. of fToll . 

a) For every g e G we have gx ~ 'Wciir], go-, g*~ l e). 

b) The Laplace transform of "Wnirj, cr, e) is 

E[e~ Hzx) ] = det(e + 2crz)^ v exp{--tr(e)}exp{^tr(e + 2cr-k zT l e)). 

c) The natural exponential family generated by "Wa(q, e, e) is {'Wo.il, cr, e) : cr e fl}. In particular, 
the statistical model 

{IVfiOu, cr, e) € P(D) : Ou, cr, 6) e / x £1 x Q} 

is an exponential family. 

d) E[jc] = rjcr + cr ★ €. 

For the remainder of this section we state and prove two propositions that we will need later in 
Section [5] The following proposition is proved in fill . 

Proposition 4.2. Suppose x ~ e\) and y ~ 'Wcj(?72 ) e, ei) are stochastically independent. 

Then the density of z = y _l ★ x is given by 



l 



• det(z)3'"-7 det(e + z )<hn^\m) eX p{_I tr ( ei + e2 )\ 



l l l l _1 1 1 

■\F\{-m + -T]2, -772, *{e + z) )oFi(-7]i, -ei ★ z)ln(z). 
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Proposition 4.3. Suppose x ~ r Wa(T]\,e,e\) and y ~ ^0(772, <?, ^2) are independently distributed. 
Then the density of u = (x + y) _1 ★ x is given by 

1 ■ det(w)5'"-7 det(e + uf^ 1 ^ exp{--tr(ei + e 2 )} 



1 l)[ + I/2 1 _i 11 

•l^iWi + 2^2, — 2 — ' 2 €l *( e + u ^ )oFi(-m, * z)ln(w)- 
Proof. If we set s = (x + y) and w = u|s, then by Corollary 14. 11 and the property a) above, we get 

s ~ < Wn(.m + m> e > 6 ) and w ~ < H / n ('7i^~ 1 )- 
The rest of the proof is similar to the proof of Proposition 14.21 □ 



5 Maximal invariants over symmetric cones 

Our primary goal in this section is to introduce a maximal invariant statistic for a family of Wishart 
distributions parameterized over an irreducible symmetric cone Q. In reference to our discussion on 
Jordan algebras we can think of Q as being the open cone of positive elements of a simple Euclidean 
Jordan algebra /. This identification, necessarily, depends on our choice of an element e e Q. which 
becomes the identity element of J. 

To restate our assumptions, let J be a simple Euclidean Jordan algebra of dimension n, rank r 
and the Pierce constant d given by the formula n = r + r(r - l)d/2. Let e be the identity element 
and Q. the corresponding open cone of positive elements of /, i.e., Q. = {x 2 \x e J, det(x) + 0}. As 
before, K is the stabilizer of e in G and c\, C2, . . . , c r is a Jordan frame in J. Recall that the group K 
acts transitively on the set of Jordan frames, therefore, for each x e J there exists a k € K such that 

X = k(gid + ■■■ +% r C r ), (5.1) 



where %\,...,i; r are the eigenvalues of x. The decomposition of x in Equation (15.11 ) is called the 
polar decomposition of x with respect to the Jordan frame c\,...,c r . The polar decomposition is 
the main tool in deriving the joint density of the ordered eigenvalues of a random vector in O. Here 
we need the following Theorem from (9l- 

Theorem 5.1. Let / be an integrable function on the Jordan algebra J. Then 

f f(x)dx - c f ( r f(k V fa) Uifj - ^) d dfiK(k))d^ ■ ■ ■ <% r , (5.2) 

where hk is the normalized Haar measure on the compact Lie group K, cq a constant depending 
only on the Jordan algebra J and R> - {(£ , . . . , £.) e W : £ > • • • > £.}. 

Notice that if x is a random vector with values in Q and var(x) + 0, then the ordered eigenvalues 
of x, almost surely, is a random vector with values in 

R; = {(£,...,&) ER r :fi >■••>£.>()}. 

Using the integral form (15.21 ) we can write a general formula for calculating the probability density 
of the ordered eigenvalues of x. Consider the proper action of K on the symmetric cone O and the 
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mapping Eig : Q — > W + (x i — > Eig(x) = (£i(x), . . . £ r (x)). Clearly, Eig(x) is a maximal invariant 
map under the action of K on Q. The distribution of Eig(x) is the joint distribution of the ordered 
eigenvalues of x. Let v denote the G-invariant measure on Q. with the density det(x)~~ln(x). By 
replacing the integrable function / in Theorem l5.1l with /(x)det(x)~~ 1q(x) and applying the integral 
formula (15.21 ) one can see that the density of the quotient measure v/fiK, with respect to the Lebesgue 
measure d£\ ■ ■ ■ d% r on R r + , is 

7=1 j<i 

On the other hand, by using Equation (7) in [2], we get the density of Eig(x) with respect to the 

quotient measure v//uk is given by Ylj=i €j J K Px(k Ey=i €j c j)dMK(k), where p x is the density of the 
random vector x. Combining these two, we obtain the density of Eig(x) as 

coY[(€j-td d f P*(kJ]ZjCj)dn K (k). 

j<i JK j=\ 

It remains to compute the constant cq, which, from Theorem 15. II does not depend on the choice of 
/. Therefore, using the Selberg integral ifTBl . 

r(jc + u - mwz + 1) 



... I exp{- ^ aj} Y\ a x f l Y\ \aj - aj\ 2z da x . . .d r = ]~] 

«-*0 :_1 ;_i ,v; :_ i 



7=1 JM 7< f 7=. 

(27r)"-T(f) r 

we obtain Co —. . In summary, 

Tn(y) 

Theorem 5.2. Suppose x is a random vector with values in Q. and p x the density of x. Then the joint 
density of the ordered eigenvalues of x , i.e., £i(x) > • • • > £ r (x) > is given by 

(2*)»-Tn$' , , . , x . . 



r Q (^) 



2 

The following Corollary is proved in H. 

Corollary 5.1. Let x be a random vector with the Wishart distribution "Wnirj, cr, e). 

a) If e = 0, i.e., the distribution is central, then the density of the maximal invariant Eig(x), thus the 
joint density of the distribution of the ordered eigenvalues of x, is given by 

Uj^j-^UU^ i , 

co i J Fo(x,-cr ). 

2sT Q (!f)det(<r)2 , » 2 

In particular, if cr = ge for some £ > 0, then this density becomes 

b) If cr = £e, then the density of the maximal invariant Eig(x) is given by 

nj<i(fj-zd d nufj h ~~ r i i ii 

co i ex Pf-^7 A €j} exp{--?r(6)} Fi(-^, x, — e). 
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Remark 5.1. Equation (15.31) is essentially proved in ifTTIl . 



Example 5.1. According to the classification of irreducible symmetric cones we have: 
i) If Q = PD W (R), the cone of N x N positive definite matrices over R, then n - N(N + l)/2, r - N 
and d = 1 Thus the density of the joint distribution of the ordered eigenvalues of the real Wishart 
distribution "Wn(t], (In), with shape parameter rj > N - 1 and scale parameter (In, is given by 



i mi * n n 

n 2 , l 



or 

n 2 ( m ' 1 



F -^exp { --g 61 n ( ^-"n*-f.). 

where r^(?7) = tt ? II/Li TC7 _ ^0 _ 1)) i s tne multivariate gamma function [ 19]. 

ii) If Q. = PDat(C), the cone of positive definite matrices over C, then n = N 2 , r = N and d = 2. 
Therefore, the density of the joint distribution of the ordered eigenvalues of the complex Wishart 
distribution CWn(tj, (In), with shape parameter rj > 2N - 2 and scale parameter (In, is given by 

, ^ ; — ex P {-i 2 ft) n f ^* - 

iii) If Q = PDat(H), the cone of positive definite matrices over H, then n = N(2N - 1), r = N and 
d - 4. Therefore, the density of the joint distribution of the ordered eigenvalues of the quaternion 
Wishart distribution WWn{t], (In), with shape parameter rj > 4N - 4 and scale parameter (In is 
given by 

si \2N 2 -2N 1 N N 1 

(2^)2^r n (i77)rn(2A0 2 tp ) = \ J y 

iv) If Q is the Lorentz cone L„, then the dimension of the Jordan algebra is n, r = 2 and d = n -2. 
Therefore, the density of the joint distribution of the ordered eigenvalues of the Lorentz Wishart 
distribution I/W„(t], (e), with shape parameter rj > n - 2 and scale parameter (e, is given by 

r(^) 1 11 



exp{-— (ft + ft)}(ftft)^"(& - ft)" 



(2^r(i77)r(^-2^)(n-3)! 2£ 

v) If n = PD3(0), then n = 27, r = 3 and J = 8. Thus the density of the joint distribution of the 
ordered eigenvalues of the octonion Wishart distribution OW(rj,(e), with shape parameter 77 > 16 
and scale parameter (e, is given by 

36(20"?'' 1 r i„ q n 1 

rri, )rx .,- 4)r( .,- 8)11!7! gftHfi^ II ■ " «?■ 

We can also use Proposition 15.11 to compute the joint density of the eigenvalues of the beta 
distribution. 
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Proposition 5.1. Let x ~ "Wo.{r]\,e, ei), y ~ 1Vn(?72, e, €2) and z = y 1 ★ x. 

a) If 61 = and €2 = e, then the probability density of Eig(z) is given by 

tinijtlu 2V2) j Ki j=x j=l L 

■iFi{^-^,\m,{e + z)-\ l -€), (5.4) 

where £1 (z) > • • • > £ r (z) are the eigenvalues of z. 

b) If £1 = e and 62 - 0, the density of the maximal invariant Eig(z) is given by 

co - -1 1 ; - n^- - ^ d fl 4'""" ft 1 + ej>~*^ ex p{-5^)} 

Bn(5J7i,5 l ») 7<! }=i >i 2 

•0^1(^1, z, \e). (5.5) 

Proof, a) Applying Theorem |5.2| to the density of the random vector z given by Proposition 14.21 the 
density of Eig(z) is then 



co- 

"~\.?>II>2'I2J j Ki j=1 /=l 

„ r m + m 1 , ^ x-i 1 x 



The proof of part b) is similar. □ 
Corollary 5.2. Suppose x ~ < Wn(^\,cr, ei) and y ~ "Wni^, o~, ^2)- 

a) If ei =0 and 62 - e, then the joint density of the distribution of the ordered eigenvalues 
£1 > • • • > Cr, of x with respect to y, is given by 

•i^i( — = — ' 2 m ' ^ e ^ '2 

&) If ei = 6 and 62 = 0, then the joint density of the distribution of the ordered eigenvalues > 
• • • > of x with respect to y, is given by 



"ifdhs n«» - w n n<> + ^ -p<-^» 

tSaijTH, 2^2) 7<; - y = i y =1 Z 



1 1 

■qFi{-t]u z, -e). 

Proof. The eigenvalues of x with respect to y are the roots of the polynomial p{£) = det(x - £y), 
which are identical with the roots of the polynomial p\{0 = det(cr _1 ★ x-£cr~ 1 *y). Therefore, we 
may assume that x ~ "Wnirji, e) and y ~ "Wn(r]2, e,cr ★ e). Now the eigenvalues of x with respect 
to y are the same as the eigenvalues of y -1 ★ x and the result follows from Proposition 15. II □ 
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The following lemma applies directly to the testing problem that we will discuss in the next 
section. 



Lemma 5.1. Suppose x ~ r Wr i (r/i,o-) and y ~ r Wn(r]2, o~, e) are two independently distributed ran- 
dom vectors with values in CI. Then the joint density of the distribution of the ordered eigenvalues, 
l\ > • • • > l r , of the random vector r = (x - y)/2 with respect to the random vector s = (x + y)/2 is 
given by 



j<i j I 7=1 



Sn( 5 77i, 2m) 

1 m + m 1 1 _i 1 

exp{--tr(e)} 1 F 1 ( — - — , -t] 2 , j(s - r) * s , -cr ★ e). 



Proof. Let £i > • • • > £ r > be the ordered eigenvalues of x with respect to y. Hence there is a g e G 
and a Jordan frame c\, • • • , c r in / such that gy - e and gx = Y!j = \ £j c j> where >•••>£,. are the 
eigenvalues of x with respect to y. Therefore, gs = £y=i 1/2(1 + Cj) c j an( l S r = T/j=\ 1/2(1 - Cj) c j- 
Without loss of generality, we may replace r with gr and s with gs. Thus we have 



■ 1 *r = V(i-^)c, 



Therefore, the ordered eigenvalues of r with respect to s are given by 

-K h = lzlL < ... < \zIl = lr< i 

and, consequently, under the transformation 

R.+ -> {(/!,-••,«: -l</ r <---</i<l} 

{&,■■■, tr) >"» (^— ^>--->T— 7"), 
1 + Cr l + 

we can compute the joint density of Zi Z r . The inverse of this map is 

{(/!,-•• ,l r ) ■ ~ 1 </,<•••< h < 1} -» 

1 - l r 1 - h 

M ( TT7;'---'iT77 ) ' 

with Jacobian 2'" rT}_i(l + f/) -2 * With this and from Corollary 15.21 we have the joint density of 
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l\ > ■ ■ ■ > l r is 



■ exp{-itr(e)} 1 F 1 (^^, X - m , (e + (s - r)" 1 * (s + r))" 1 , ^ * 6) 



2 



2 r 1 (771 + 772 1 1 _, 1 
2 exp{--tr(e)} 1 F 1 , -rj 2 , -(s - r) ★ s , -cr* e) 



n 7 r =1 d + /,) 



''l + ''2 r r 

BClWl'lW j<i j=\ j=l 

■ exp{-^tr(e)hFi(^-^, -772, -{s - r) * s~ l , ^cr * e), 
where we have used the fact that z = (s - r) _1 ★ (s + r). 



6 A generalization of Bartlett's test 

In classical multivariate analysis Bartlett's test is used for testing equality of variances of two, or 
more, normally distributed populations based on their observed samples. Since the sample covari- 
ance matrix of a normally distributed population follows a Wishart distribution, Bartlett's test is 
indeed a testing problem concerning the scale parameters of two, or more, Wishart models. In this 
section we use the results obtained in previous sections to show how we can set up a similar testing 
problem for scale parameter of a Wishart model parameterized over a symmetric cone. 

Suppose j] > (r - \)d is a known parameter, and consider the statistical model 
CWnfa. 0-1) ® 0-2) € P(Q x Q) : (crj , cr 2 ) eflxQ) 

and its submodel 

CWnfa. <r) ® Wnir], cr) e P(Q x Q) : a e Q). 
We wish to test the null hypothesis Hq : cr\ = o"2 = cr vs. // : <j\ + cr 2 . 

Theorem 6.1. For the observation (x,y) e Q.xQ. the maximum likelihood estimator of cr under Ho 
is s(x,y) = (x + y)/2, and the likelihood ratio statistic for testing Hq vs. H is riy=i(l _ if)* 11 ' where 
^1 > ^2 > • • • > h are the eigenvalues of r with respect to s. Furthermore, under the null hypothesis 
Ho, the statistics s(x, y) and n{x, y) = • • • , l r ) are independently distributed, six, y) ~ 'Wo.iri, cr) 
and the density of n(x, y) is given by 

Proof. Consider the diagonal action of G on Q. x Q. defined by (g, (x,y)) h-> (gx,gy). Evidently, 
jt is surjective, and invariant under G. Even more, we show that n is maximal invariant under G. 
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Suppose n(x,y) = n(x',y'). We may assume, without loss of generality, that x = x' - e, and 
therefore, the eigenvalues of y and y' are identical. Let YIj=\ lj c j an d Zy=i Ijc'j be, respectively, the 
spectral decompositions of y and y' . Since K acts transitively on the set of Jordan frames of Q, we 
have y and y' are in the same orbit of G. It is known that under the null hypothesis Hq the statistic 
s(x,y) = (x + y)/2 is the maximum likelihood estimator (MLE) of cr, which is also invariant under 
the action of G. Therefore, the independency of s(x, y) and n(x, y) follows from [3 , Lemma3]. The 
density of tt(x, y) is given by Lemma I5TT1 for rj\ = 772 = r] and e = 0. To complete the proof let 
calculate the likelihood ratio test LR. Since (x,y) is the MLE of (en, cr 2 ) we have 



q(x,y) 



Lo(c?|(x,y)) 



det(^Z)-'' 



H(o-\,a-2)\(x,y)) Aet(x)-^ d&t{yY 



2 r,r 



det^^'detCy)^ 
det(x + y) 1 ! 



det^" 1 *y)ii 
det(e + x~ x -k y)n 



r 



I _ (l+^) 2 j 



4-n 



1-// \^ 

( — -) 

v i + 1/ 



J=i (1 + T^) 2 



□ 
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